load('recalibrated_data.mat')

global delta_TS beta_TS rho_TS THETA_TS

beta_TS = 0.5030;
delta_TS = 0.9880;
rho_TS = 1.1051;


for THETA = -0.25
    THETA_TS = THETA;
    
    batchmsm20160515dr_exp_singleRun;
    load('simulation_output/hs_xz_exp_rho0_dr_singleRunsim')
    suffix = strcat('outside_avgZ_TC');
    mydata.(suffix) = avgZ_;
    suffix = strcat('outside_avgX_TC');
    mydata.(suffix) = avgX_;
    suffix = strcat('outside_avgTC_TC');
    mydata.(suffix) = avgTC_;
    
end


save('recalibrated_data','mydata')
